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CN : 1 Preliminaries 
o 

' Let B (0, i?o) C denote a three-dimensional spherical droplet of radius Rq > 0, centered at the 

■ origin. Let Sq denote the set of symmetric, traceless 3x3 matrices i.e. 



So = {Qe M3x3; Q. . = Q^.,; Q-- = O} (1) 

where M^^^ is the set of 3 x 3 matrices. The corresponding matrix norm is defined to be 



< 

\Q\^ = QijQij i,J = l...3 (2) 

and we will use the Einstein summation convention throughout the paper. 

We work with the Landau-de Gennes theory for nematic liquid crystals [7] whereby a liquid crystal 
_ configuration is described by a macroscopic order parameter, known as the Q-tensor order parameter. 

■ Mathematically, the Landau-de Gennes Q-tensor order parameter is a symmetric, traceless 3x3 

matrix belonging to the space 5*0 in ([1]). The liquid crystal energy is given by the Landau-de Gennes 
I energy functional and the associated energy density is a nonlinear function of Q and its spatial 

derivatives [3 [T7]. We work with the simplest form of the Landau-de Gennes energy functional 
that allows for a first-order nematic-isotropic phase transition and spatial inhomogeneities as shown 
below [I2] - 

Ilg[Q]= / ^^|VQ|2 + /b(Q) dV. (3) 

Jb{0,Ro) ^ 

^ \ Here, L > is a small material-dependent elastic constant, |VQp = Qij^kQij,k ( note that Qij^k = 



For our purposes, we take fs to be a quartic polynomial in the Q-tensor invariants as shown below 



' ^ '■' ) with i,j,k = 1 ... 3 is an elastic energy density and : 5*0 — )• M is the bulk energy density. 



/b(Q) = "^^/*^ rQ^ - |trQ3 + ^ (trQ^)^ (4) 

where trQ^ = QijQjpQpi with i,j,p = 1 . . . 3, a, 6^, > are material-dependent constants, T is 
the absolute temperature and T* is a characteristic temperature below which the isotropic phase 
Q = loses its stability. We work in the low-temperature regime with T « T* and hence, we can 
re-write @ as 

/b(Q) = -ytrQ2 - ^trQ3 + ^ (trQ^)^ (5) 
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where > is a temperature-dependent parameter and we will subsequently investigate the — )■ oo 
limit, known as the low-temperature limit. One can readily verify that Jb is bounded from below 
and attains its minimum on the set of Q-tensors given by jlH [T^ 



Qmin = {Qe So;Q = s+{n0n - - ] , ne S'} (6) 



where I is the 3x3 identity matrix and 



62 + + 24a2c2 
= ^2 • (7) 

We are interested in characterizing global minimizers of the Landau-de Gennes energy functional 
in ([3]) , on spherical droplets with homeotropic or radial anchoring conditions [13] • The global Landau- 
de Gennes minimizers correspond to physically observable liquid crystal configurations and hence, 
are of both mathematical and practical importance. We take our admissible Q-tensors to belong to 
the space 

A={Qe W^'^{B{0,Ro);So);Q = Qb on dB{0,Ro)} (8) 

where VF^'2 (5(0, Rq); Sq) is the Soboblev space of square-integrable Q-tensors with square- integrable 
first derivatives [5]. The Dirichlet boundary condition Q;, is given by 



Qfe(x) = S+ ( - ^ ) G Qn,in (9) 



where x E is the position vector and |^ is the unit- vector in the radial direction. For completeness, 

we recaU that the T^^'2-norm is given by ||Q||y[/i.2 = (^Jb^oRq) IQI^ + I^QP dV^ ^ and the L°°- 
norm is defined to be ||Q||/^oo = ess supxgs(o,i?o)IQ(-'^)l 13 • 

We define a modified Landau-de Gennes energy functional as shown below - 

Ilg [Q] = / ^|VQ|2 + fB (Q) - min /b (Q) dV (10) 
Jb{0,Ro) ^ Qe-^o 

where (Q) — minQgSg fs (Q) > for all Q G Sq. The existence of a global minimizer for the 
modified functional Ilg in the admissible space A is immediate from the direct methods in the 
calculus of variations [5]; the details are omitted for brevity. It is clear that Q* G ^ is a minimizer 
of Ilg if and only if Q* is a minimizer of Ilg in ^ and hence, it suffices to study minimizers of 
the modified functional in ()10p . In what follows, we study the modified functional in (jlOp and drop 
the hat for brevity. 

For a fixed o? > 0, let Q" € A denote a corresponding Landau-de Gennes global minimizer. The 
Euler-Lagrange equations associated with Ilg are given by a nonlinear elliptic system of coupled 
partial differential equations: 

LAQij = -a'^Qij - ^2 ^QipQpj - ^QpqQpgSij^ + (trQ2) Q^j i,j,p,q = 1...3 (11) 

where ^QpqQpqSij is a Lagrange multiplier accounting for the tracelessness constraint [12]. It follows 
from standard arguments in elliptic regularity that Q" is smooth and real analytic on B{0,Rq). 
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The notion of a limiting harmonic map was first introduced in [12] and is crucial in what follows. 
A limiting harmonic map Q'^ G ^ is defined to be 

Q° = «+(n°®nO-I) (12) 

where n*^ is a minimizer of the Dirichlet energy [TS] 

/[n] = / |Vnp dV (13) 

in the admissible space = |n G VF^'^ [B{0, Rq); S'^) ; n = on dB{0, Rq)^ . In the case of a 

spherical droplet with homeotropic boundary conditions, n*^ is unique and given by the radial unit- 
vector [1^. Hence, the limiting harmonic map is unique for our model problem and is given by 

We note that Q*^ has a single isolated point defect at the origin. 

In what follows, we keep L,6^ and fixed in (|10p . Following the methods in jl4j . we first 
introduce a re-scaling relevant to the low-temperature limit: — )■ oo. Let t denote a dimensionless 
temperature 

27a2c2 , , 

t = —^>Q (15) 

so that the — >• oo limit corresponds to the t — )• oo limit and define 



h+ = ~W- t^oo. (16) 



Define 



1 / 27c^ 



and the corresponding Landau-de Gennes energy functional is given by (up to a multiplicative 
constant) 



iLG[Q] = / ^IVQI^ + ^ f (1 - - ^A/|trQ3 + 



dV (18) 



'B{0,i?o) 

where L = > is fixed, f{t) is a function that can be explicitly computed, 

(l - |Q|2)^ - §^y/|trQ3 + f{t) > for Q G 5o (see definition of Ilg m i^) and for t sufficiently 
large, we have 

for positive constants cri,cJ2 and constants 71,72 independent of t in the t — )• 00 limit. The corre- 
sponding admissible space is 



A=\QeW'^^iBiO,Ro);So);Q = ^^^QbondBiO,Ro)\ (20) 
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where Qt has been defined above. In what fohows, we drop the bars from the re-scaled variables 
for brevity and all subsequent statements in this section are to be understood in terms of the scaled 
variable in (fT7|) . 

Next, we quote some results from [12] that are relevant to the development of our mathematical 
framework. 



Proposition 1 For each t > 0, let Q* = s* (n* (g) n* — ^l) denote a uniaxial global minimizer o/Ilg 
in l[T8\}. in the admissible space A defined in where s* G W^''^ {Q,R) and n* G W^''^ (0,5^). 

Then there exists a sequence {t^} with t^— T-ooasA:— )-oo such that Q**-' — )■ strongly in 
W^''^ {B{0, Rq); So) , where Q,^ is the re-scaled limiting harmonic map defined in (T^. 



Proof: The proof of Proposition [T] is identical to the proof of Lemma 3 in jl2| and the details 
are omitted for brevity. From the global energy minimality, we have 



Ilg[Q*1 < Ilg[Q°] = E{L,Ro) yk 



(21) 



where > is a positive constant independent of tk i.e. the Q**-' 's have bounded energy in the limit 
— )• oo. Note that 

t 



:i-|QT)'-^/|tr(Q°f + /(t) 



t 







for all t > 0, by definition ( since G Qmin)* One immediate consequence of the strong convergence 
result are the following integral equalities - 



B{0,Ro) 



|VQ°p dV = lim 



lim 



:i-iQ*N) 



B(0,iffl) 

2 8h. 



From the inequalities ()19p . we deduce that 



t 



VQ**|2 dV 
y|tr(Q*'=)^ + /(t) 



dV = 0. 



[ (1- |Q*^|2)' = 0. 

JBiO.Rn) 



lim 



(22) 



(23) 



□ 



Proposition 2 For each t > 0, let Q* denote a global minimizer of Ilg 'ii^ U8\) . in the admissible 
space A defined in [2^) . Let {t^'} be a sequence such that t^ ^ oo as fc —>• oo. Then 



lim IQ^Hx)! < 1 Vx G B{0,Ro) 

k—^oo 



(24) 



Proof: The proof follows from a standard maximum principle; see [12] Proposition 3 for an analogous 
statement with proof. □ 



Proposition 3 [I12j; Lemma 2] For each t > 0, let Q* denote a global minimizer of 1i,g in il8\) . 
in the admissible space A defined in Ii20\) . Define B{x,r) = {y G B{0,Rq); |x — y| < r} C B{0,Rq) 
and 



e(Q*,VQ*)=^^|VQf + ^ 



:i-iQ 



t|2\2 



8/li 



tr{Q'f + f{t) 
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Then 

- [ e(Q*,VQ*) dV<^ [ e(Q*,VQ*) dV, VxGS(0,i?o), r < R, (25) 

^ JB{x,r) ^ Jb{x,R) 

SO that S(x,i2) C B{0,Ro). 

Proof: The proof can be found in [[l2]; Lemma 2]. An analogous boundary monotonicity formula 
can be found in Lemma 9 [12]. □ 



Proposition 4 For each t > 0, let Q* denote a global uniaxial minimizer of Ilg in- 1^18]) . in the 
admissible space A defined in Let {t''} be a sequence such that t oo as A; oo with 

Q*fe — ). Q'' in W^''^ {B{0, Rq), So) as k ^ oo, where Q,^ is the re-scaled limiting harmonic map 
defined in [T4\ )- For any compact K C B{0,Ro) such that K does not contain any singularities of 
Q*^ i.e. does not contain the origin, we have 

lim |Q*Hx)| = 1 Vx G (26) 

and the limit is uniform on K. 

Proof: Proposition m is a consequence of the pointwise uniform convergence 



lim 



:i-|Q*N)'-^/|tr(Q*^)' + /(t.) 



everywhere away from the singular set of Q*^ i.e. away from the origin. This uniform convergence 
result holds in the interior and up to the boundary. The proof can be found in [T^, Propositions 4 
and 6. □ 

To summarize, let {t^^ be a sequence such that — t- oo as /c — t- oo and let {Q*''} denote a 
corresponding sequence of purely uniaxial Landau-de Gennes minimizers. Then up to a subsequence, 
Q*fe -f QO in W^'"^ (-8(0, Rq), Sq) as A; — )■ oo. Further for k sufficiently large, IlgIQ**"] can be bounded 
independently of tj- since 

IlgIQ*^ <Ilg[Q°] 

and IQ*''! is strictly positive (and bounded from below) everywhere away from the origin. 

2 Statement of main results 

Our main result is the following : 

Theorem 1 Let B{0,Rq) C M'^ denote a spherical droplet of radius, Rq, centered at the origin. For 
each > 0, let Q" denote a global Landau-de Gennes minimizer in the admissible space A defined 
in 13^. Then in the limit a} — )• oo, Q° cannot be purely uniaxial i.e. cannot be of the form 

Q"(x) = s(x) ( n®n- ^ ) x£B{0,Ro) (27) 



3, 

for a function s : B{0, Rq) — t- M and a unit-vector field n G W'^''^ (-^(0) -^o)) S"^) ■ 
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Comment: From Jl]/, it is known that i/ Q = s (n (X" n — ^) E W'^''^ {Q,So) on a simply- connected 
domain Q C R^, then n G W^'"^ (O, 5^) . 

We prove Theorem [T] by contradiction in the subsequent sections. Of key importance is the 
division trick used in [18] for the Ginzburg-Landau theory for superconductivity in three dimensions. 
We adapt the division trick in |18j to the Landau-de Gennes framework for nematic hquid crystals. 
We point out the following important differences: (i) we have a parameter and we are interested in 
the asymptotics of global energy minimizers in the — )• cxd hmit; the Ginzburg-Landau equations in 
|18j are parameter- free, (ii) the nonlinearities in the Landau-de Gennes equations (jlip are different to 
the nonlinearities in the Ginzburg-Landau equations introducing additional technical complexities 
and (iii) the Landau-de Gennes macroscopic variable is a two-tensor field Q G Sq whereas the 
Ginzburg-Landau macroscopic variable in M.^ is a three-dimensional vector field u S M^. 



To make better contact with the framework used in [18] , we introduce the following dimensionless 
variables as in 



where t = ^'^^4'^^ > is the reduced temperature [12], t > 1 throughout the paper and /i+ has been 

defined in (|16|) . The length-scale = We note that the position vector x has been re-scaled 

in (|'28|) whereas it was left unchanged in (|17p . The corresponding dimensionless energy density is 



e(Q,VQ) = -|VQ|^--trQ^-^^trQ^ + ^(trQ^) +C{t) (29) 

h 

where C(t) = 2 + ^ — w additive constant that ensures 



--trQ^-^^trQ^ + ^(^trQ^j +C{t)>0 VQ G 5o. 
The corresponding Landau-de Gennes energy functional is given by 



Ilg[Q] = / e(Q,VQ) dV, (30) 

JB{0,Rt) 

where Rt = 7\/ti?o for a fixed constant 7 > 0. In particular, i?t —?• 00 as t —)• 00. In what follows, 
we drop the tilde on the dimensionless variables for brevity and all subsequent results are to be 
understood in terms of the dimensionless variables. Prom ([8]) and ()20p . the admissible Q-tensors 
belong to the space 



(31) 



= |q G (5(0, Rt),So) ; Q = yi - ^l) on dB{0, i?i)| 

The associated Euler-Lagrange equations are [131 H] - 

AQij = -Qi, - (^Q^kQkj - Ytr(Q2)j ^ ^Q^.tr(Q2), ij = 1,2, 3. (32) 

All global and local energy minimizers in are classical solutions of (j32p . 
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Proposition 5 For each t > 0, assume that a uniaxial global minimizer of T^g exists in the ad- 
missible space Aq defined in we denote this uniaxial minimizer by Q*. Let [t^] be a sequence 
such that t'^ —7- cxo as A; —7- oo with Q*'' — )• Q*^ in W^'"^ {B{0, Rt^.), Sq) as k ^ oo, where is the 
limiting harmonic map 

Then from Proposition [7| 

lim / Atr^^-^tr^^ + ^(tr^A\cit,)=0 (33) 

*fe^°° yB{0,Rt J ^ tk ^tk V / 

where Rt^ oc y/tk- For k sufficiently large, we have 



(i) Q*fc = -^flQ'H (n®n - il) for some n G W^'^iB{0, Rt,), S'^) . 

(a) Q*'' is a classical solution of the following nonlinear system of elliptic partial differential equations 
on -6(0, Rt^) where Rt^ oo as tk ^ oo : 



AQ*) = (|Q*N' - 1) (^^ + ^ {W - |Q*N) Q| (34) 



(Hi) |Q*'=|(x) < 1 for allxe B{0,Rt^), 

(iv) ;^Ilg[Q*'=] < 127r and 

(v) Ym^tk^^ Q^KO) = 0. 

(vi) All derivatives of Q^'' can be bounded independently of tk in the /c — t- oo limit i.e. 

l|V^Q*MlL^(B(0A,))<C^i i>i (35) 
for a positive constant Cj independent of tk- 

Proof: In what follows, we drop the subscript k for brevity and work in the t — )• oo limit. Proof 
of (i): This follows from the uniaxial character of Q* i.e. it has two equal eigenvalues 

Q*(x) = A(x) (e(x) e(x) + f(x) ® f (x)) - 2An(x) ® n(x) 

where e, f , n form an orthonormal frame at x and n(8)n + e®e + f(8)f = 1. Using the above, Q* 
can be written in the simpler form 

Q*(x) = -3A(x) (^n(x) ® n(x) - ^I^ (36) 

where 

IQT(x) = 6A2(x). (37) 
From [11], a global uniaxial Landau-de Gennes minimizer has A < and from ()36p -(|37p. 
— 3A = Y^IQ*I- The representation formula in (i) follows. □ 

Proof of (a): If Q* is a uniaxial global Landau-de Gennes minimizer, then it is a classical solution 
of ()32p . The partial differential equations ()34p follow from substituting the representation formula 
in (i) into ([32]). □ 

Proof of (Hi): The proof follows from multiplying both sides of by Q*^ and applying a standard 
maximum principle argument for |Q*P; the details are omitted for brevity. □ 
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Proof of (iv): This is a direct consequence of the global energy minimality. The limiting harmonic 
map Q'' is simply given by 



in terms of the dimensionless variables in ()28p . A direct computation shows that ^^^•[Q'^] = 127ri?j 

(since — ^t: 
and hence 



(since -ifr (Q^)^ - ^^^tr {Q^f + ^ (tr (Q°)^)^ + C{t) = by definition because Q° € Q 



mm J 



Xlg[Q*] < Xlg[Q°] = l^TTRt. (39) 

The inequality in (iv) follows. □ 

Proof of (v): This follows from Proposition^ We have a topologically non-trivial boundary condition 
Qfe in ([9]) and hence every interior extension of must have interior discontinuities. The extension 
n in (i) has interior discontinuities and at every such point of discontinuity x*, Q*(x*) = (see [H] 
for further discussion on these lines; Q* is analytic at x* whereas n is not and n can lose regularity 
only when the number of distinct eigenvalues of Q* changes. The number of distinct eigenvalues of 
Q* can change only when Q* relaxes into the isotropic phase i.e. Q*(x*) = 0.) From Proposition HI 
as t — 7- oo, all isotropic points are concentrated near the singular set of Q*^ and the singular set of 
qO j-Qgj-gjy consists of the origin. Hence, we have limt^ooQ*(0) = 0. □ 

Lemma 1 Assume that Q* is a uniaxial global Landau-de Gennes minimizer in the admissible space 
Aq on the droplet B{0, Rt) in the limit t — )■ oo. Then the bulk energy density satisfies the following 
inequality - 

Proof: From Proposition [5l we have < |Q*| < 1 and one can check from the representation formula 
in Proposition [5] (i) that 

Substitute the above into the definition of /(Q*) in (P0|) 



/(Q') = 4lQf->f + |lQr4 + ^-| (41) 



and one can check from ()4ip that 

□ 

Our second main result concerns the characterization of uniaxial global Landau-de Gennes min- 
imizers if they exist. 

Theorem 2 Assume that Q* G is a uniaxial global Landau-de Gennes minimizer on the droplet 
B{0, Rt) in the limit t — )• oo, where i?t oc i/t — oo as t ^ oo. Then Q* is an entire solution of 

AQ*,. = (IQf - 1) Q*,. + ^ (IQf - |Q*|) Q*, (42) 
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with Q*(0) = and limf_>.oo ■^^^^[Q*] < 12-k. There exists a T G 0(3) such that 



Q*(x) = ^(|x| 



Tx(g)Tx 1 

3^ 



(43) 



where h : [0, oo) — t- M"*" is the unique, monotonically increasing solution of the ordinary differential 
equation 



d'^h 2 dh 6h 
dr^ r dr 



3h. 



-h+ ^ 
t 



{h^ - h^) 



(with r = subject to the boundary conditions 

h{0) = h{r) ^1 r ^ oo 
< h{r) < 1 

in the limit t — )■ oo. 



(44) 



(45) 



The last ingredient of the proof of Theorem [T] is the foUowing result from [151 1 
Theorem 3 1151 \14^ : The radial-hedgehog solution 

H = /i(|x|) 



X (X) X 



il 

3 



(46) 



where h is the unique solution of |^^[ ) subject to |^5[ j, is a solution of the Euler- Lagrange equations 
[34\ )- Moreover, H is not a global Landau-de Gennes minimizer in the admissible space in the 
limit t ^ oo. 



Proof: One can explicitly construct a biaxial perturbation of the form 

H.(x) = H(x) + ^;:^(i-^) (.«.-i) 

with 0" = 10 and z = (0,0, 1), that has lower Landau-de Gennes energy than the radial-hedgehog 
solution in the limit t — )■ oo, Rt — )• oo. This perturbation does not have an isotropic core at the 
origin. The details can be found in |14t I15j. □ 

To summarize, we obtain symmetry results for global uniaxial Landau-de Gennes minimizers 
in the limit t — )• oo i.e. we show that every uniaxial Landau-de Gennes minimizer is the radial- 
hedgehog solution in (|46|) (modulo a rotation) in the limit t — )• oo. The final step is to use results on 
the radial- hedgehog solution in the limit t — t- oo from |15tll4j to prove Theorem [TJ In the subsequent 
sections, we proceed with the proof of Theorem [2j 



3 Symmetry of uniaxial Landau-de Gennes minimizers 

The main trick is to adapt the division trick of [18] to the Landau-de Gennes framework for nematic 
liquid crystals. Let Q* denote a global uniaxial Landau-de Gennes minimizer in the space Aq for t 
sufficiently large so that Proposition [5] holds. Define 

S«(x) = |W (47) 

where h is the unique solution of (j44|) subject to (|45p in the t — >• oo limit. We first prove some 
auxiliary properties of the function h : [0, oo) — t- [0, 1) in ()44p . 
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Theorem 4 Define the radial-hedgehog solution H as in J^fip ; H is a solution of the system of 
partial differential equations for all t > 1. In the limit t — >• oo, 
(i) h : [0,00) — [0, 1) is analytic; 



Hij(O) = and VHij(O) = 0, 
(ii) < h{r) < 1 and 



h{0) = h{0) =0;h (0) > 0; 



^2 ^2 

< h{r) < 



r2 + 14 



r2 + 



t ' 



(Hi) we have the following gradient bound 



< Ci||VH||^c»™3) < C; 



M|x|) 



where Ci and C2 are positive constants independent oft as t ^ 00, 
(iv) we have the following hound for the second derivative 



h (|x|) < C3||V2H||i^(K3) < C4 



Hi, 



where C3 and C4 are positive constants independent of t as t — )■ 00 and 
(v) we have the following hound for the third derivative (noting that 



1; 



h (|x|) 



h 



< C5IIV H||2,cx)(K3) < Cq 



(48) 



(49) 



(50) 



(51) 



(52) 



where C5 and Cg are positive constants independent oft as t —t- 00. 

Proof of (i): The analyticity of /i, the relations (|48|) and Hjj(O) = have been proven in 
To prove VHjj(O) = 0, we use the following equality 

Hi,-(x)Hi,-(x) = /i2(|x|), i,j,= 1...3 (53) 

so that for any fixed direction characterized by the unit- vector e^ where a = 1 ... 3, we have 

H,,(x)H,,-,(x) = /i(|x|)/i'(|xr''" 



where x„ = x • e^. We set x = |x|eci and re- write the above as 

lim H.,(|x|e ) -H.,.(0) \kiM) - 

|X|^.0 |X| J. M I l^l^g 



M|x|) 



(54) 



which implies Hjj^Q(0)Hjj^Q,(0) = /i'(0) = for each a = 1 ... 3. Hence, VHjj(O) = 0. □ 
Proof of (ii): The proof of (ii) can be found in [15] and [14j . 

Proof of (Hi): We start with the relation (j53p and differentiate both sides with respect to the 
coordinate direction e^ to find 



h{\A) 



Hjj Hjj Q, Xq, 



(55) 
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Here 





< 1 and 




|x| 




h 



1. From ()55p . we have that 



^ Ci||VH||^oo 



M|x|) 



(56) 



where Ci > is independent of t. Prom the gradient bound derived in Lemma A.l in [2] and 
standard results in elhptic regularity (also refer to Proposition [5] (vi)), we have that ||VH||^cx)(ir3) 
and all higher derivatives of the radial-hedgehog solution H can be bounded independently of t in 
the 4—7-00 limit. The inequality (j50p now follows. □ 

Proof of (iv): We start with the relation (|53p and differentiate twice to obtain 



dh dh , d^h 
+ h- 



5xq Ox^ d:x.ad:x.i3 
A direct computation (see shows that Hjj^cHjj^^ = -§^-§^ so that 



(57) 



h 



(58) 



and hence 



< CallV^Hllio 



< Ca 



for positive constants C3 and C4 independent of t in the i — )• 00 limit. Pinally, it suffices to note 
that 



M|x|) 



and (|5ip now follows. □ 

Proof of (v): We compute an explicit expression for q^^^qI^q^ as shown below 



h (|x|): 



M|x|) 



X^(5a^ + Xa(5, 



/37 



2xq,x^x^ 



X 



+ 



+h (|x|): 



|x| 



XaX/3 



|X 

We multiply both sides of ([59]) by 



h (|x|) 



+ /i(|x| 
- to obtain 



Sap^-f + ^a-f'^p + <J7/3Xq, 3Xq,X/3X, 



+ 



(59) 



XofX^X^ 



9xa3x^9x^ 



Straightforward computations show that 



9Xq,(?X/35x^ 

and on combining (|60p and (|6ip. we obtain 



/i (|x|) 



HjjHjj^Q,^-y Xq,X^X^ 



(60) 



(61) 



(62) 



The inequality (|52p now follows. □ 
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Lemma 2 Assume that Q* is a global uniaxial Landau-de Gennes minimizer in the admissible space 
Ad defined in (EIP with Q*(0) = 0. Then VQ*(0) = 0. 

Proof: From the results in jllj . we have that a uniaxial Landau-de Gennes minimizer can be written 
in the form 

Q* = y|lQ*l(n^5n-il 

for some n G VF^'^(i?(0, S*^). Therefore, |Q*| has a minimum at the origin and we have 
V|Q*|(0) = 0. Using the relation Q-^Q-j = |Q*P , Q*(0) = and V|Q*|(0) = 0, we can repeat 
the same steps as in Theorems (i) (see ^ and dMI) to deduce that VQ*(0) = 0. □ 

Proposition 6 Assume that Q* is a global uniaxial Landau-de Gennes minimizer in the space Aq^ 
for t sufficiently large, so that Proposition holds. Define 

/i(|x|) 

Then 

|S(x)|<C7 VxG5(0,i?j) (64) 
where Cj is a positive constant independent of t. We have the following bounds for VS ; 

f |x|2 < 1 

|VS(x)|< W ' ' - (65) 

[Cio |x|2 > 1 

where Cg and Ciq are positive constants independent oft. 



Proof: From Theorem U we have 



so that 



lx|2 

|X| , 2 



/i(|x|)>^ |x|^<l (66) 



h 



< i^lQ 



X|2 



where ||Q*||loo(m3) < 1, Q*(0) = and VQ*(0) = (see Lemma[2|). From standard Taylor expansion 
formulae [?], we have 



Q*^.(x) = Q*^. „^(sx)x„x^(l - s) ds |x|2 < 1 

^0 



(67) 



so that 

|Q*,(x)| < ||v2q*||^oo(k3)|x|2. (68) 

Q* is a classical solution of ()34p and from [2] and standard results in elliptic regularity (also 
see Proposition [5] (vi)), we have that ||VQ*||loo, ||V^Q*||loo and all higher derivatives of Q* can be 
bounded independently of t in the f — )• oo limit. Substituting the bounds ()66p and (|68p into the 
definition of S, we obtain 

|S(x)| < Cii |x|2 < 1 (69) 

where Cn is a positive constant independent of t. For |xp > 1, it suffices to note that ||Q*||_loo(ir3) < 1 
and /i(|x|) > ^ for |x|2 > L Hence, |S(x)| < for |x|2 > 1 where C12 is a positive constant 
independent of t in the t — t- 00 limit. The inequality ()64p now follows. 
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Next, we compute bounds for VS. Consider the case |xp < 1 first. We re-write S as 

. QV|x|2 



so that 



IVSI < D 



+ 



M|x|)/|x|2 
iQVlxpl 



|V(Q7|xf 

/i(|x|)/|x|2 ' (/i(|x|)/|x|2)2 



|V(M|x|)/| 



X 



(70) 



(71) 



where D > is a constant independent of t. We start with the integral formula (|67p and compute 
V (Q*/|xp). Straightforward but tedious computations show that 



V 



< -A X p < 1 



X 



where Di > is independent of t. Similarly, we have for |x|2 < 1 (see Theorem U]), 



M|x|) 



so that 



V 



M|x| 



h {s) (|x| — s) ds 



D2 , ,9 

< t4 |x| < 1 

X 



(72) 



(73) 



(74) 



where D2 > is independent of t. From the bounds in Theorem 2] and ()68p . we have that and 
can be bounded independently of t for jxp < 1. Substituting ([72|) and ([7^ into ([7T]) . we have 



|VS(x)| < ^ |x|2 < 1. 



For |xp > 1, we have 



dS 



9x^ 



h 



x^ 
|x| 



and ||/i'(|x|)||ioo(iR3) can 



For |xp > 1, /i is bounded away from zero, ||Q*||l°°(ir3) ^ 1; II^Q*IIl° 
be bounded independently of t. Hence 

|VS| < L>4 |xp > 1 

for a constant D4 > independent of t. The bounds ()65p now follow. □ 

We need more careful estimates for S and its gradient near the origin for the subsequent analysis. 



Proposition 7 Let Q* denote a global uniaxial Landau-de Gennes minimizer in the space Aq^ for 
t sufficiently large, so that Proposition\^ holds. Define S as in \6^) . Then 



Xq-X^ 

«J ~ "ijaP I |2 



Sjj(x) = B, 

d 



5x^ 



9x^ 



B, 



- + o(l) |x|^0 

XqX/j" 



X 



+ 0(1) |x| 



(75) 



where 'Bija/^ 



/i"(0) '^'-^'"'/^ 



1 ... 3 is a constant matrix. 
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Proof: The proof of Proposition [7] follows from lengthy and technical computations. We state 
the main points for completeness. An explicit computation shows that 



dSjj _ _d_ 



dy 



^tjaP I |2 



where 



A*^.^ (Q*; h) + S*.^ (Q*; h) + C*-^ (Q*; h) 



%/\^? d //i(|x|) 



(76) 



(M|X|)/|X|2)2 5X^ V |xP 



|xp d 

/i(|x|) ^y^J 
d 



Qii,a/3(0) 



XgX^ 

2|x|2 



5? 



Qji,a/3(0) 



X^vX 



2|x|2 



M|x|) /."(O) 



(77) 



The next step is to show that A* ,B* , C* are bounded as |x| — ?• in the t ^ oo limit. 
We recall that 



for |x|2 < 1 (see Theorem S]) and 



Q*(x) 



X 



h 1 

^ - 15 



where is a constant independent of t (see (j68p ). A straightforward computation shows that 



/i(|x|) = /i"(0)^+ /'"'r(r). 



|x|-r)^ 



dT 



(78) 



so that 



d fh{\x\)\ _ d //i(|x|)\ x^ 



dx^ \ |xp 
X. 



(9|x| \ |x 



X 



T \ T 



u7 Mr) 1 - rr 02 



X 



The equality (|79|) implies that 



and from (|52|) . we have 



d fhi\x\) 



d: 



x^ 



< DeWh \\lo 



d fhi\x\) 



dxj \ |xr 



< D7 



(79) 



(80) 



(81) 



where Dq, Dj are positive constants independent of t. The inequality ()8ip immediately implies that 

|^:,^(Q*;/i)| <i?8 (82) 



where Dq > is independent of t. 
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Next, we turn to B*j^ (Q*;/i). Recalling that Q-j(O) = VQ\j{0) = 0, we have 



21x1 



2x 



Therefore, 
M|x|) 



9 / Xq-X/jx^ 



2|x|2 



9x^ V 2|xp 



gli?;,,(Q*;/^)= / (1-.) 
x| Jo 

Hence, 

(Q*;/i)| < Dg [||V4Q*||i.o(K3)|x| + 1] < Z^io 

where Dg and Dio are positive constants independent of t. 

Finally, we turn to C*j^ (Q*;/i). From ()78p . we obtain the following inequality 



(83) 

ds. (84) 
(85) 



x|2/i"(0)-2/i(|x|) 
which when combined with the inequalities 

d ( 



< Dn\\h 



(86) 



9x^ \ |x| 



< Rll and /i(|x|) > ^ for IxP < 1 
X 15 



yields 

(Q*;/i)| <^i3 

where Dii,Di2 and Z^is are positive constants independent of t0. 
From ([82]), ([85]) and dSTD, we deduce that 



(87) 



dSij _ d 



B, 



XgX^ 



0(1) as |x| —7- and t — ?■ oo. 



Next, we consider an explicit Taylor expansion for Sij near the origin: 



b,,ixj - |x|2 +^^^"^2/i(|x|)|x|^ 



|x|^/i"(0)-2/i(|x|) +^y" (l-s)2Q*.„^^(sx)x„x^x^ds. 

(89) 



From (1861). we obtain 



qX^ 



2/i(|x|)|x|^ 



/i (0) -2/i(|x|) 



< Dl4||V2Q*||icx,m3)|x| < L>i5|x 



and 



1 /"^ 



< Z)i6||V^Q*||ico(iR3)|x| < Di7|x| 

where D14 — Diy are positive constants independent of t. Hence, 

Xq,X/3 



<Di8|x| |x^0 



(90) 



for a positive constant Z^ig independent of t. Proposition [7] now follows. □ 



^From Theorem U h (0) is a constant independent of t for t sufficiently large. 
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Theorem 5 Assume that there exists a global uniaxial Landau-de Gennes minimizer Q* G for 
each t > 0. Let {t^} be a sequence such that — )• oo as A; — )■ oo so that (up to a subsequence), 
— ^ Q*^ in W^'^ {B{0, Rt^), So) as tj. — )• oo, where Rk oc and is the unique limiting 

harmonic map defined in i38\). Define 



(Skh (x) 



Q-)(x) 
M|x|) 



i,j = 1...3 



(91) 



as in i63\) . where h : [0, oo) — )• [0, 1) has been defined in ^4\ ) ' i45\ )- Then in the limit k — )■ oo, 
SkeC^ {R^\{0};So), 



lim — - 

k^oo Rf^ 



lim -— 



1 



1 



fc^oo Rk JB{0,Rk) |x| 



4 

dV = 0. 



< 127r 



(92) 



Proof: In what follows, we drop the subscript k for brevity and work with t^ large enough so 
that Propositions [5l [6] and [7] hold. From global energy minimality and the inequality established in 
Lemma [H we have (see Proposition [5|) 



1 

'Rt 



1 



VQf + IQf ^y^ 



t JB{0,Rt) 



1 

R't 



1, 



vqT + /(Q') < 12^ 



t JB{0,Bt) 



(93) 



where i2t oc -v/t and /(Q*) has been defined in Lemma [TJ From the strong convergence in Proposi- 
tion [1] and the inequality established in Lemma [H we have 



lim — 

t^oo Rt Jb 



{o,Rt) 



We note that 



|2|2 



< 2 



t^oo i2i JB{0,Rt) 
t|2|2 



|1 - IQ* 

h^ 



(94) 



(95) 



For |xp < 1, |S| < Gi for a positive constant Gi independent of t; see Proposition [6l For |xp > 1, 



we recall the bounds in Theorem S] to find 

2 



|1 - 



< 



Si 



X 



|2|2 



< <52 1 



|x| > 1 

-IQfl' 



1x1^ > 1 



where 61,62 > are independent of t. Combining the inequalities (|95|) and (j96j) . we obtain 



lim — 

t-^oo Rf^ J^, 



1 



< 6:^ lim 



1 



' i-s>oo R- 



[ (1-IQf)' dV = 

■t JB{0,Rt) 



(96) 



(97) 



KO,Rt) 

where > is independent of t. 

Next, we turn to the elastic term, |VSp. For |xp < rg where tq >> 1 is a constant independent 
of t, we use the estimates in Proposition [6] to obtain 

1, 



S(0,ro) 



|VSr dV < 6ir^ 



(98) 
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where (54 > is independent of t. On the region B{0,Rt) \ i?(0, tq), an expHcit computation shows 
that 



|VS|2 = ^^ + |Qf 



/l2 



For |x| > ro, we recah from [T3] that 

h'{\^\) 

for a positive constant ^5 independent of t and 

ivsp 



< 



VQf (1 + -^) 



(99) 



(100) 



X 



for a positive constant 6q > independent of t, in this regime. Combining ([97j) - (jlOO|) . we deduce 
the following chain of inequalities 



lim 



1 



1, 



ivsr + 



1 - isi 



dV = lim 



< lim 

t--s>oo Ri 



■t JB{0,Rt) 
1 



1 



t-^oo R^ JB{0,Rt) 2 



VSP dV < 



|VQ*p < 127r 



B{0,Rt) 



(101) 



where the last inequality follows from Proposition [SJ 

Finally, we turn to the integral lim(_>.oo J^^q j^^-^ ^^"|^p ^ dF. Recall that |S| satisfies the global 
upper bound |S| < 67 for a positive constant 67 independent of t (see Proposition [6]) . Then a direct 
computation shows that 

/ ' I ' ' ' dV < 5s (102) 

for a positive constant independent of t. On the region B{0, Rt)\B{0, 1), we use Young's inequality 
to deduce 



|1 - ISI 



1 

< - 
- 2 



;i-is|2)^ 4 



X ■ 



(103) 



Combining ()102p - (|103p and recalling ([97]) . we obtain 



lim 



1 



t-s>oo Rf 



|1 - IS| 



dV < 



t JB(0,Rt) 



X 



< lim 

t— >oo 



+ 



1 



2\2 



:i-isn 



dV + 



Rt '^Rt J B{Q,Rt)\B{0,l) 4 2Rt 

as required. The proof of Theorem [5] is now complete. □ 



dV 



B{0,Rt)\B{0,l) 



X 



(104) 



3.1 PDE-methods for S 

Let Q* G denote a uniaxial global Landau-de Gennes minimizer for t sufficiently large so that 
PropositionsElEllTland Theorem [5] hold. Define S as in ([63]) . We recall that Q* is a classical solution 
of 

3h. 



AQ 



IQ 



t|2 



1 + 



t 



(IQf -|Q*I) 



Q 



(105) 
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and /i is a solution of the ordinary differential equation 



h" + l-h' --^ = h'-h+^ (h' - h') 
X xr t ' 



(106) 



where /i' denotes etc. Straightforward but lengthy manipulations show that S satisfies the 
following system of partial differential equations - 

AS., + (l + ^] (1 - S., = -2^S.,,,g-^ + ^/.S,, (1 - |S|) i, j = 1 . . . 3. (107) 



Following the methods in [18], we multiply both sides of the partial differential equation ()107p 

(108) 



with Sjj^fcj^- One can check that 



1 /5S 



+ 



5 



9x„ 



x| " |x| \(?|x 
1 + ^) /i'(|x|)(l-|S|2)S,,S,,-fc 



2 X X 



1 + 



2\2 



(1-|S|^) 



2/i/i + 



^ „ Xi, x„ 
O __G1 i_ 9 



h f 



|x| 

2^ 

r 

2 



9 /x„/i2(l-|SP)'^ 



9Xp \ |x| 



9 



x|2 ^-''^Ixl 5x„ 



/i V5|x| 

Sl(i-isp] 



and 



(109) 

(110) 
(111) 



Using (fTUSj) - (fTTTD . we obtain 



5 



9> 



2- ■ ixi 



^ 3/i+\ Xp /t^(l-|Sp)^ 3xp(l-|S| 

t y |x| 4 |xp 



1 /SS,; 



x| V^|x| 



+ 1 + 



3/i+\ (1 - |S 



2^2 



2hh + 



2/i2 



+ 2- 



/i \d\^\ 



3h. 



■/i(l- |S|)S,,Si#^) 



We recall from Theorem U] and the results in |14j that h' > for |x| > 0, for t sufficiently large. 
Define 



|x| 



/ 3/i+\ x„/i2(l-|S|2)" 3x„(l-|Sn ^ ^ 



Lemma 3 VFe /iaue 



^p^dA ^ 127r as (5 ^ 



'|x|=5 |X| 

where dA is the surface area element on the sphere of radius 5 centered at the origin. 
Proof: By the definition of <I>p in ()113p . we have 

<^r,^dA = 



(114) 



|x|=<5 |X| 



|x|=<5 



-|VS|^ 



dS 

d\ 



x|y V * / 4 |x| 



(115) 
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From the estimates in Proposition [71 we have the fohowing as 5 — )■ 

2 



IVSI 



V 



dS 



1 - IS 



o(|x| 



2 l-'^l^ |BjjQ,/3XQ,X^| 



+ 0(1) 



(116) 



where Bj^Q,^ 



is a constant matrix and i,j,a,(3 = 1 ... 3. 



Substituting ()116p into ()115p . we get 



|x|=5 |x| 



1 



|x|=5 



V 



X 



+ 3 



— |Bjjo/3XQ,X^|^ 



+ 0(1x1-2) 



127r + o(l) + 



|x|=5 



V 



X 



A direct computation shows that 



1 



|x|=<5 



V 



X 



X 



BijQ,/3XcjX/3 



= 



for the constant matrix 'Bija/s 
Lemma 4 The integral 



h"iO) 



and hence, the conclusion of Lemma [3] follows. □ 



1 



|x|=l 



V 



BjjQ/jXQ-X^ 

W2 



Bjjci/3XQ,X^ 



X 



dA = 



for any constant Bij^p such that Bijajs = Bjiajs, Bij^ji = Bj^-^a and Bij^a = Bu^p = 0. 



(117) 



(118) 



Proof: Consider the matrix 



B 



tja/z 



h"{0) 



i,j,a,/3 = 1...3. 



(119) 



Recalling that Q* is a symmetric, traceless 3x3 matrix which is a classical solution the PDE-system 
(jlOSp (so that we have Qjj = Qij sa ^om equality of mixed partial derivatives), we obtain 



Bijap — Bjjfj^, BjjQ^ — Bjj^Q,, BjjfjQ, — 0, BjjQ^ — 

for all i,j,a,(3 = 1 ... 3. 

A direct computation shows that 



(120) 



V 



BjjQ^XQ,x^ 



(121) 
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and 



so that 



J|x|=l ^ 



V 



X 



X 



dA 



26, 



rp 



dA 

|x|=l J|x| = l 



for i,j,p,q,r,s = 1 ... 3. 

Using spherical coordinate representation, we can check that 



|x|=l 



■Jqs 



and 



I XpXqX^XgdA — — — [Spq6rs ~l~ ^pr^qs ~l~ ^ps^qr] • 

J|x|=l 



Substituting (fT2i|) and (fT25]) into (fT23|) . we obtain 



|x| = l 

47r 
"3 



V 



-Bjja^XQ,X^ 



' 3 




|x|2 


|x|2 



dA 



r, [2^ijrs'^ijrs ^ijpp^ijss ^ijqr^ijrg '^ijsr'^ijrs 



(122) 



(123) 



(124) 
(125) 



(126) 



and the right-hand side vanishes by virtue of the properties estabhshed in ()120p . The integral 
equality (jllSp now follows. □ 

Proposition 8 Let Q* € denote a uniaxial global Landau-de Gennes minimizer for t sufficiently 
large so that Propositions {^{^^ and Theorem\^hold. Define S as in h6S^) . Then 



xeB{0,Rt 



|S(x)| = l xGB{0,Rt) 
where Rt = /^V^ and fi is a constant independent oft. 



(127) 
(128) 



Proof: We integrate both sides of (|112p from |x| = to |x| = Rt, divide by Rt, use Lemma [3] 
and take limit t — )• oo to obtain 



1 



Rt 



127r + lim — , , ,,.011 
t^oo Rt Jq Jb{o,r) |x| V^I^I 



1 / dS 



+ 1 



3/ii 



(i-isn 



r 2/^2- 
2hh + — 



+ 



2h ( as 



h \d\-s\ 



dVdR 



— lim 



3h. 1 



t^oo t Rt 

: lim — / 1|VS|2 

t^-oo Rt JB{0,Rt) 2 



/ ' / /i(l-|S|)S,jS,j-fc^ dV dR = 
Jo Jb{o,r) m 

SV / 3h+\h^l-\S\^)' 3(1-|S|2) 



dS 

d\2 



dV. 
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From (19211 ■ we have that 



1 f 3(1-|S|2) 



hm ^ / ' ' dV = 



and 



+ (1 + ^1 — ^ — < 127r 



since /i^(|x|) < 1 on B{0,Rt) and hm(_j.oo = 0- 

From Theorem m we recall that h is monotonically increasing in the limit t — t- oo so that every 
term in the integrand of 



i_ f^^ r 1 /dsy / 3h^\ (i-|sp)^ 



'B{0,R) 
is non-negative. 



2hh + 



2h^ 



+ dVdR 



h V5|x| 



Finally, we estimate the integral J^^ Ib{o R) ~ l^l)^«i^«i>fc^ ^ follows: 



B{0,R) 

< C 



[ h{l - \S\)S,,S,,,k^ dV 
Jb(o.r) 1^1 

/ {i-\s\rdvVl[ 

Jb(o.r) Jb 



IVSP dV 



IB(0,R) 



1/2 



<C*./f[R\VR (130) 



where lim/j_>>oo "^-^ = 0. for all i? > where C and C* are positive constants independent of R. The 
first inequality follows from Cauchy-Schwarz and the fact that |S| < D for a positive constant D 
independent of R (see Proposition [6]) . From Theorem [SJ we have 



lim — 

-R— !>oo it 



/ (1-|S|2)' (iF = 



so that if we set 



/ (1-|S|2)' dV = f[R], 



then f[R] = o(i?) as — t- cxo. The second inequality also uses the upper bound Jbi^or) |VSp < 
127ri2 from yielding the inequality (|13Up . 

Substituting the upper bound (|13U|) into (|129|) . we get 



hm 

t— >oo t 



^ [ ' [ HI- |S|)SySy-fc^ dV dR 
Jo Jb{o,r) m 



(131) 



since -^Z f[R]R dR = o{R^) as i? — oo and ~ as t — > oo. 



Combining the above, we deduce that 

/ (i-isiy 



2/i /i 



2/i2i 2/i' / as \^ 



X 



+ 



h \d\^\ 



dV = Q (132) 



for all ii > and the conclusion of Proposition [8] now follows. □ 
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Proposition 9 Let Q* G denote a uniaxial global Landau-de Gennes minimizer for t sufficiently 
large so that Propositions El El and Theorem\^ hold. Then in the limit t — )• oo, we have 

Q* = /i(|x|) - X G B{0, Rt) (133) 

where h is the unique solution of (44\ ) subject to (J^ ) and Rt = ii\ft, with fi being a constant 
independent oft. 

Proof: From Proposition [51 we have that in the Umit t — )• oo 



X 



X 



Q* = M|x|)M,, — (134) 



where Mjj = y f (™ ^ m — |) for some m G (from the uniaxial character of Q*), 

|M(x)|^ = l xG5(0,i?f) (135) 
and Mjj — )• (^xf ~ ^) ^ ^'^^ ~^ ^ from the imposed Dirichlet boundary condition. 
We substitute ()134|) into the governing system of PDEs for Q* as shown below - 



AQ 



IQf -1) + ^(|Q*P-|Q*|) 



Q\j, (136) 



multiply both sides of (|136p with Mjj (noting that MjjMjj ^ = 0) to find 



r 

Consider the minimization problem 

|2 



VM|2 = 3|Vm|2 = — . (137) 



min / |Vn|^ dV (138) 



B{0,R) 



where B{0, R) is a three-dimensional droplet of arbitrary radius R> and 



X 



iV = <^ n G W'^^{B{0, R); S^); n = — on dB{0, R) } . (139) 



X 



It is known that the minimization problem (|138p - (|139|) has a unique minimizer [TU] 



nmin = 1^ (140) 

and 

/ iVriminl^ dV = SttR. (141) 

Jb{o,r) 

The unit-vector field m G W^^'^{B{0, Rt); 5^), m = |fy on dB{0, Rt) and from (fT37)l . we have 

|Vmp dV = SnRt. (142) 



'B{0,R) 

Comparing ()14ip and (|142p . we deduce that m is a minimizer of the problem ()138p - (|139p on -6(0, Rt] 
and from the uniqueness of the minimizer, we deduce that 

m(x) = 1^ x£B{0,Rt). 
The conclusion of Proposition [9] now follows. □ 
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